arXiv:hep-ph/9808369v2 16 Oct 1998 


Radiative decays of D mesons 

S. Fajfer Q , S. Prelovsek a and P. Singer 6 

“ J. Stefan Institute, Jamova 39, P.O. Box 300, 1001 Ljubljana, Slovenia 
b Department of Physics, Technion—Israel Institute of Technology, Haifa 32000, Israel 


The short-distance contribution due to the c —> wy penguin as well as long-distance contributions are 
considered in an analysis of D —> Vy decays, the latter being dominant. The matrix elements for nine 
D —> Uy transitions are calculated using a hybrid model which combines heavy quark effective theory 
and the chiral Lagrangian. We present the expected range of the branching ratios for these decays, the 
most frequent ones D° —> A*°y and D+ —> p + y, being a few times 1(D 4 . 


1. INTRODUCTION AND THE Q <?y 
TRANSITION 

Radiative flavour changing quark transitions 
Q —> q 7) occur at the one loop level in the elec- 
troweak Standard Model and would show up at 
four different scales s —> dy, c —> ity, b —> s(d )y 
and t —> c(u)y. These transitions are of intensive 
theoretical and experimental interest, especially 
since it has been point out [1] that the QCD en¬ 
hancement of b —> sy brings it into the realm 
of observability [2], Here, we are concerned with 
the c —> uy transition, which can be expected to 
induce hadronic decays like D —> Vy. 

The amplitude for the transition of a 

heavy quark Q to a light quark q and an on-shell 
photon is given by [3] 
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to an inclusive branching ratio of 10~ 8 only [5]; 
accordingly, one expects for the exclusive transi¬ 
tions D —> Uy a branching ratio of about 10~ 9 
from this short-distance contribution. 

The long-distance (LD) contribution to c —> ity 
has also been estimated recently [6]. Consider¬ 
ing the process c —> u + fifi with the quark 
pairs fifi hadronizing into vector mesons, one de¬ 
rives the amplitudes for the c —> uVi processes 
(Vi = p,u>,<p) using the QCD corrected nonlep- 
tonic weak Lagrangian. Taking then the trans¬ 
verse part of the c —> uVi amplitude, with gauge 
invariance and vector meson dominance one ar¬ 
rives at the amplitude [6] 
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where 02 is a Wilson coefficient and Cv MD is given 
by 


where in our case F 2 = Y,\V\ C V£ U F 2 : \ and the 
summation is over A = d,s, b. F 2 is calculable 
in the electroweak model [3]; furthermore, the 
inclusion of QCD corrections leads to a strong 
enhancement of its value [4,5]. Despite this in¬ 
crease, the short distance c —> tty amplitude leads 
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Using the vector-photon couplings measured in 
leptonic decays, and assuming gv t { 0) ~ gvA^v ) 
one finds [7] a strong cancellation in (3) due to 









the GIM mechanism and SU(3) symmetry, lead¬ 
ing again to small branching ratios, e.g. about 
10 ~ 8 for D° —> p °7 [ 6 ], only slightly larger than 
the SD contribution. 

2. A MODEL FOR D —► V 7 
TRANSITIONS 

The results of the previous section indicate the 
need for a detailed treatment of the LD contribu¬ 
tions to the D —> V 7 transitions; indeed, several 
approaches were put forward using a pole model 
[4], a quark model [ 8 ] and QCD sum rules [9]. 
In the present treatment, we aim for a more com¬ 
prehensive and systematic treatment for these de¬ 
cays, employing an effective hybrid Lagrangian 

[ 10 ] which combines two approximate symmetries 
of QCD, the infinite heavy quark (Q = c ) mass 
limit and the chiral limit for light quarks. 

We assume [6,11] the radiative decays to orig¬ 
inate from the nonleptonic weak transition D —> 
VVq [4,12] followed by the conversion Vq —■> 7 
via vector meson dominance. In addition, there 
are transitions due to direct photon emission from 
the initial D state. The effective nonleptonic La¬ 
grangian we use is given by 
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and in order to evaluate the matrix elements of 
the product of two currents we use systematically 

[11] the factorization approximation. The quark 
bilinears in (4) are treated as interpolating fields 
for the appropriate mesons, the relevant hadronic 
degrees of freedom being the charm pseudoscalar 
( D ) and vector mesons (D*) and the light pseu¬ 
doscalar (P) and vector mesons (V). In the fac¬ 
torization approach we use, the D —> VV 0 ampli¬ 
tude is schematically approximated as 
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There are three terms in the factorized D —* VVq 
amplitude, referring to an annihilation part, a 
Vo-spectator part and a U-spectator part, each 


of them involving several diagrams. For the 
QCD-induced constants 01,02 we take ai = 1.26, 
02 = —0.55 as determined from an extensive ap¬ 
plication of (4) to nonleptonic D decays [13]. 

The general gauge invariant amplitude for the 
decay D(p) —> V(pv)"/{k) is 
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We have classified all diagrams contributing to 
Ape, Apv, and the explicit expressions are given 
in references [ 6 ] and [ 11 ], 

The calculation requires the {V\ (V ^ — A M ) D) 
matrix element, which has one vector V{q 2 ) and 
four axial-vector Ai(q 2 ) form factors. Actually, 
gauge invariance and finiteness reduce the num¬ 
ber of independent form factors to three. For 
these we take a pole dominance behaviour, with 
known values of pole masses. Information on 
E(0), A\ (0), A 2 (0) we extract from the semilep- 
tonic decays D + —> K*°£p e and D\ —» ip£v e . 


3. RESULTS 


The calculation of the decay amplitudes in¬ 
volves various constants of the hybrid La¬ 
grangian. Their imprecise knowledge is the first 
source of uncertainty. In addition, the ampli¬ 
tudes contain several components with unknown 
relative phases, which is yet another source of 
uncertainty. This allows us to give only an ex¬ 
pected range, for each of the calculated transi¬ 
tions. In the Table below, we give the predicted 
range of the calculated values for the branching 
ratios. The first two decays are Cabibbo-allowed, 
the next five are Cabibbo-forbidden and the last 
two are doubly forbidden. To give an indication, 
the photon energy in the first two decays is 717 
and 834 MeV respectively. As the table shows, a 
few of these decays are in a range which should 
become experimentally feasible soon. Their de¬ 
tection will contribute decisively to the under¬ 
standing of long-distance dynamics in D decays. 

On the other hand, we are still left with the 



question of how to detect directly the c —> wy 
transition. One possibility is related to acciden¬ 
tal cancellations between the various other contri¬ 
butions to the amplitudes of Cabibbo-forbidden 
decays [ 6 ], which is probably unlikely. The alter¬ 
native would be to search for weak radiative tran¬ 
sitions among baryonic states containing several 
c quarks, like S++ -> £++7 and fi++ -> S++ 7 . 
In these cases, due to the valence structure of the 
participating baryons, the c —> uj process would 
play a dominant role[14]. 

Table 1 

The predicted branching ratios for nine D —> V 7 
decays 


D 

■ P 7 Transition 

Br Ratio x 10 5 [11] 

D° - 

-4 A'*°7 

6-36 

Df 

-» P + 7 

20-80 

D° - 

-»■ P° 7 

0 . 1-1 

D° - 


0.1-0.9 

D° - 


0.4-1.9 

D+ - 

-»■ p + 7 

0.4-6.3 

Df 

-4 K *+ 7 

1.2-5.1 

D+ - 

-4 K * +r ) 

0.03-0.44 

D° - 

-4 K *° 7 

0.03-0.2 
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